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We determine, by numerical simulations on a lattice, the gauge-invariant two-point correlation functions of 
the gauge field strengths in the QCD vacuum with four flavours of dynamical staggered fermions. 
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The gauge-invariant two-point correlators of 
the gauge field strengths in the QCD vacuum are 
defined as 



r (x) = <0|Tr {G w (x)SG ury (0)S^} |0) , (1) 



where G w = gT a Gp\ p is the field-strength tensor 
and S = S(x, 0) is the Schwinger phase operator 
needed to parallel-transport the tensor G va (0) to 
the point x. 

These correlators govern the effect of the gluon 
condensate on the level splittings in the spectrum 
of heavy QQ bound states • They are the ba- 
sic quantities in models of stochastic confinement 
of colour and in the description of high- 

energy hadron scattering (see Ref. Q] and refer- 
ences therein). 

A numerical determination of the correlators 
on the lattice already exists in the quenched (i.e., 
pure-gauge) 5/7(2) theor y m , and also in the 
quenched 5/7(3) theory |p|,|l0|. Here we present 
results obtained in full QCD, i.e., by including the 
effects of dynamical fermions JL1[ . Four flavours 
of staggered fermions and the Wilson action for 
the pure-gauge sector have been used. The de- 
termination has been done at two different values 
of the quark mass. 

In the Euclidean theory, the most general 
parametrization of the correlators is the follow- 
ing §§: 

T^IJ,p,vu{x) = (SftuSpcT - S^aSpv) [D(x 2 ) + T>i(x 2 )] 
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where T> and T>\ are invariant functions of x 2 . It is 
also convenient to define the following quantities: 
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On the lattice we can define an operator T>^ ip va , 
which is proportional to V^p^a- in the naive con- 
tinuum limit, i.e., when the lattice spacing a 

a 4 V l± {d 2 a 2 ) + 0(a 6 ) . 
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However, the naive continuum limit of Eq. (Q) is 
spoiled by the presence of lattice artefacts, i.e., 
renormalization effects from lattice to continuum 
due to the short-range fluctuations at the scale of 
the UV cutoff. In order to remove these artefacts 
we adopt the same technique used in Refs. P,|l0|]. 
The basic idea is to remove the effects of short- 
range fluctuations on large distance correlators by 
a local cooling procedure 0Jl3|] . 

We have measured the correlations on a 16 3 x 24 
lattice at distances d ranging from 3 to 8 lattice 
spacings and at [3 = 5.35 ((3 = 6/g 2 , where g is 
the coupling constant) . We have used a standard 
Hybrid Monte Carlo algorithm, in particular the 
so-called ^-algorithm described in Ref. (w) . The 
bare quark mass was chosen to be a ■ m q = 0.01. 
A determination was also made at a ■ m q = 0.02. 

The scale of our system is fixed by the physical 
value of the lattice spacing a. We shall use the 
following parametrization 

231/625 

exp 



-25^1.(5) 



where the scaling function /(/3) = Ap ■ a((5) is 
given by the usual two-loop expression for gauge 
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Figure 1. The functions D_l/A^ (upper curve) 
and T)\\/A F (lower curve) versus physical dis- 
tance, for quark mass a ■ m q = 0.01. 
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Figure 2. The functions T>±/A F (upper curve) 
and T>\\/Ap (lower curve) versus physical dis- 
tance, for quark mass a ■ m q — 0.02. 



group SU(3) and Nf — 4 flavours of quarks. Ap 
in Eq. (H) is an effective A-parameter for QCD in 
the lattice renormalization scheme, with Nf = 4 
flavours of quarks. At the value of (3 used in our 
simulation the lattice spacing, extracted from the 
string tension or the p mass, isa~0.11±0.01fm 
|l5| , [il| , so that our lattice size is approximately 2 
fm and therefore safe from infrared artefacts. 
In Fig. 1 we display the results for V»/A F and 



T>±/Ap versus d 



phys 



d ■ a, for a quark mass 



a ■ m q = 0.01. As in Rcf. [Q, we have tried a 



best fit to these data with the functions 



V{x z ) = A, exp - 
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The continuum lines in Fig. 1 correspond to the 
central values of this best fit (see Ref. [p|). 

The corresponding results for the quark mass 
a ■ m q = 0.02 are displayed in Fig. 2. From 
a ■ m q = 0.01 to a ■ m q = 0.02 the effective Ap 
does not change appreciably within the errors, so 
that we have assumed the same value of Ap as for 
a ■ m q = 0.01. Again, the continuum lines in Fig. 
2 correspond to the central values of the best fit 



with the functions (||) (see Ref. |iT[). 

A quantity of physical interest which can be 
extracted from our lattice determination is the 
correlation length Xa of the gluon field strengths, 
defined in Eq. (^|): it is relevant for the descrip- 
tion of vacuum models . To obtain the value 
of Xa in physical units, the physical value of the 
lattice spacing must be used. This gives, for the 
first quark mass a ■ m q = 0.01: 

X A = 0.34 ± 0.02 ± 0.03 fm {a-m q = 0.01) ,(7) 

where the first error comes from our determina- 
tion and the second from the uncertainty in fixing 
the physical scale. 

Similarly, for the second quark mass a ■ m q = 
0.02 we obtain: 

X A = 0.29 ± 0.01 ± 0.03 fm {a-m q = 0.02) .(8) 

The values (Q) and (0) must be compared with 
the quenched value ]9|,|10| 

X A = 0.22 ± 0.01 ± 0.02 fm (YM theory) . (9) 

The correlation length Xa decreases by increas- 
ing the quark mass, when going from chiral to 
quenched QCD. Of course, a determination of the 
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physical value of Xa should be done with more 
realistic values for the quark masses. 

Another quantity of physical interest which can 
be extracted from our results is the the so-called 
gluon condensate, defined as 

G 2 = : G; U G; V :) (a s = . (10) 

As first pointed out by Shifman, Vainshtein and 
Zakharov fT?j] , it is a fundamental quantity in 
QCD, in the context of the sum rules. 

The gluon condensate can be expressed, in 
terms of the parameters defined in Eq. (^), as 
follows El: 



6 

G 2 ~ -^{Aq + A^ 
At a ■ m a 



(11) 



0.01 this gives, in physical units, 



Go = 0.015 ±0.003 



+0.006 
-0.003 



GeV 4 



(12) 



At a ■ 



0.02 we obtain: 



Go = 0.031 ± 0.005 



+0.012 



GeV 



-0.007 VJCV ■ (13) 

These values should be compared with the corre- 
sponding quenched value JTo| : 

G 2 = 0.14±0.02 GeV 4 (YM theory) .(14) 

As expected, the gluon condensate G 2 appears 
to increase with the quark mass, tending towards 
the (pure-gauge) value of Eq. (14). We can try 
to understand the dependence of G 2 on the quark 
masses using the following low-energy theorem 
|l8| , valid for small quark masses: 

j I a s . pa pa ,\ _ ^4 

where b = 11 — |^V/, for a gauge group SU(3) 
and Nf quark flavours. For a ■ m a = 0.01 we have 
approximately m/ ~ 44 MeV Making use 

of the popular values for the quark condensate 
((qq) ~ -0.013 GeV 3 H,@) and for the phys- 
ical quark masses (m u ~ 4 MeV, nrid — 7 MeV 
and m s — 150 MeV), we can extrapolate from 
the value ( |l2"| ) to the physical gluon condensate, 
obtaining the following estimate: 



: ( : If If 



(15) 



G 



(physical) 



0.022 GeV 



(16) 



The prediction ( |l6| ) agrees with phenomenologi- 



cal determinations fll9|,|2C 
0.011 GeV 4 . 



G 



(empiric) 



0.024 ± 
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